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ABSTRACT 

In 1992 the Far-Ultraviolet Space Telescope (FAUST) provided measnrements 
of the nltraviolet (140-180nm) diffnse sky backgronnd at high, medinm, and 
low Galactic latitndes. A significant fraction of the detected radiation was 
fonnd to be of Galactic origin, resnlting from scattering by dnst in the diffnse 
interstellar medinm. To simnlate the radiative transfer in the Galaxy, we 
employed a Monte Garlo model which ntilized a realistic, non-isotropic radiation 
field based on the measnred flnxes (at 156nm) and positions of 58,000 TD-1 
stars, and a clond strnctnre for the interstellar medinm. The comparison of 
the model predictions with the observations led to a separation of the Galactic 
scattered radiation from an approximately constant backgronnd, attribnted to 
airglow and extragalactic radiation, and to a well constrained determination of 
the dnst scattering properties. The derived dnst albedo a = 0.45 ± 0.05 is 
snbstantially lower than albedos derived for dnst in dense reflection nebnlae and 
star-forming regions, while the phase fnnction asymmetry g = 0.68 ± 0.10 
is indicative of a strongly forward directed phase fnnction. We show the 
highly non-isotropic phase fnnction to be responsible, in conjnnction with the 
non-isotropic UV radiation field, for the wide range of observed correlations 
between the diffnsely scattered Galactic radiation and the colnmn densities of 
nentral atomic hydrogen. The low dnst albedo is attribnted to a size distribntion 
of grains in the diffnse medinm with average sizes smaller than those in dense 
reflection nebnlae. 


Subject headings: radiative transfer — scattering — nltraviolet: ISM 
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1. Introduction 


When illuminated by the galactic interstellar radiation field, dust in the interstellar 
medium gives rise to scattered radiation known as diffuse galactic light (DGL). Efforts to 
observe the DGL at different wavelengths and attempts to derive scattering properties of 
the dust grains from such data have extended over the past half-century, and corresponding 
work directed at the far-ultraviolet spectral region has been ongoing during the past 


quarter-century. A detailed review of DGL studies was given by Witt (|1990|) , while Bowyer 


( p.991|) and Henry (|1991|) reviewed in-depth the particular complications involved in the 
observation and interpretation of the far-ultraviolet background radiation due to several 
sources, including the DGL. 

Major challenges facing DGL studies in the far-ultraviolet include the need for 
observations with high diffuse-source sensitivity which provide extensive sky coverage while 
simultaneously offering the means to separate the flux from discrete sources, such as stars 
and galaxies, from that of the diffuse background. For the interpretation of DGL data, there 
has been a demand ([Mathis 1993|) for multiple-scattering models which take into account 
the inhomogeneity of the scattering interstellar medium and the anisotropic distribution of 
the interstellar radiation field, which is extreme for the far-ultraviolet sky (jMurthy fc Henry 
1995[). The successful flight of the Far-Ultraviolet Space Telescope (FAUST) ([Bowyer et al 


1993[ ), conducted as part of NASA’s 1992 ATLAS-1 shuttle mission ([Graven fc dorr 1988[ ), 


has been a response to the first of these two challenges. The photon-counting imaging 
detector of FAUST ([Lampton et al. 1990[) allowed a separation of the diffuse background 
from that due to stars in the field to a level where unresolved stars contribute less than one 
percent of the total diffuse intensity in the field ( [Cohen et al. 199^) . At the same time, this 
instrument provided sensitive measurements of the diffuse background in fields covering 
over 1000 square degrees in representative directions ranging from high to intermediate 
latitudes ([Sasseen et al. 1995[) . Sasseen et al. (1995) showed that the dominant component 
in the measured far-ultraviolet background is due to starlight scattered by galactic dust by 
demonstrating a clear relationship between spatial power spectra of the IRAS lOO^um cirrus 
and the FAUST far-ultraviolet diffuse background images. 

An initial attempt to model the FAUST ultraviolet background and to constrain 
the scattering properties of the diffusely distributed dust was undertaken by Sasseen 
& Deharveng ([1996[) . This attempt, which used the simple model of Jura ( [1979[) for 
high-latitude dust illuminated by a constant plane source, met with only limited success 
however. The model was unable to correctly reflect the anisotropy and the longitudinal 
variation of the intensity of the illuminating radiation field. 


In this paper, we are using a substantially more sophisticated model to treat the 
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radiative transfer for the diffuse background regions observed by FAUST. This model 
was hrst introduced and described briefly by Witt & Petersohn (|1994| , hereafter WP) in 
the context of an analysis of ultraviolet background observations made by the Dynamics 
Explorer 1 spacecraft ( [Fix et al. 1989 ). 

The organization of this paper is as follows. In §2 we will summarize the data reduction 
and present an initial analysis. This will establish that the measured intensities are a 
combination of DGL, which is dependent on the column density of the dust in the line of 
sight, and of roughly constant contributions from airglow and possible extragalactic diffuse 
background radiation. In §3 we will present the WP model in greater detail then was done 
in WP (|1994|) . Following this, we carry out the radiative transfer analysis of the FAUST 
data in §4, followed by a discussion of the implications of the derived dust properties for 
the likely processing of interstellar dust in the diffuse medium away from the galactic plane 
in 55. In 56 we will summarize our results and formulate our conclusions. 


2. The Data 


The data used for this analysis were obtained with the imaging telescope, FAUST, 
during a March 1992 shuttle flight. FAUST has a bandpass covering the range 140 - 
ISOnm and a 7.6° circular held of view. The detector was a photon-counting micro cha nn el 
plate with a wedge-and-strip position sensitive anode ([Lampton et al. 1986|) . The in-hight 
angular resolution of the all-rehecting camera (D=161mm, f/1.12) was approximately 3.8', 
which permitted the identihcation and removal of point sources down to a detection limit 
of 1 X 10“^^ erg s“^ cm“^ A“^. Each pixel found to contain excess hux due to a point 
source was replaced with a weighted average based on distance, of those surrounding pixels 
containing negligible hux from the source (Sasseen et al. 1995). 


A summary of the data used for this analysis is given in Table 1. The list contains 
the same FAUST helds as studied by Sasseen & Deharveng (|1996|) , to which we have added 
data from the FAUST image ^ 2 in Dorado. The columns in Table 1 give, respectively, (1) 
the name of the FAUST held (Sasseen Sz Deharveng 1996); (2) the FAUST image number 
( jBowyer et al. 19^31) ; (3) the number of 0.5° x 0.5° pixels in each image for which the dihuse 
UV intensity was determined; (4) & (5) the galactic coordinates for the image center, (6) 
the average intensity of the dihuse radiation in the 140-180nm FAUST band in units of 
photons cm“^ s“^ sr“^, (hereafter referred to as ’’units”); (7) the average intensity 

of the dihuse IRAS background radiation at 100//m ([IPAC 1991|) ; and (8) the average HI 
column density taken from the Bell Lab HI Survey ( Stark et al. 19921) , supplemented with 
data from the Parkes HI survey (jCIeary et al. 1979| ; pdeiles fc Cleary 1979|) for declination 

























Table 1. 


Data Used in this Study 


FAUST Field 

# 

Points 

1 

h 

<I(UV)> 

<I(100/ini)> 

<N(HI)> 




(°) 

(°) 

(PU) 

(MJy/sr) 

(lO^^/cm^) 

(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

(8) 

Dorado 

2 

49 

266.2 

-44.4 

1705± 90 

0.272±0.227 

200.8± 55.1 

NGPKM 

8 

64 

250.1 

76.3 

900± 56 

2.055±0.435 

261.5± 17.5 

M87 

9 

64 

285.5 

75.0 

1431± 75 

2.389±0.367 

242.8± 17.8 

Virgo P2 

10 

58 

284.3 

76.4 

1023± 51 

2.255±0.445 

245.5± 18.0 

Virgo PI 

11 

62 

283.8 

71.6 

1052T137 

2.315±0.241 

191.6± 29.8 

Centaurus 

14 

59 

302.8 

21.7 

3169±174 

8.624±1.104 

747.1T 76.0 

M83 

16 

64 

314.7 

32.8 

1777T147 

4.878±0.536 

446.0± 52.4 

NGC 6752 

18 

62 

336.8 

-25.2 

2190±289 

4.263±0.862 

539.0± 50.3 

Hydra 20a 

20 

62 

261.4 

37.6 

1376±221 

3.470±0.531 

539.8± 69.3 

Hydra 20b 

20 

49 

267.5 

38.1 

1211T184 

3.054±0.619 

478.5± 37.1 

Hydra 20c 

20 

64 

273.5 

38.5 

1361±228 

2.550±0.610 

450.2± 36.6 

Hydra 21a 

21 

63 

314.5 

35.3 

2200±227 

5.409±0.673 

539.0± 63.2 

Hydra 21b 

21 

50 

309.5 

38.3 

2379±241 

8.165±1.209 

774.7± 72.5 

Hydra 21c 

21 

62 

304.3 

40.7 

2179±210 

7.063±1.152 

709.7±103.3 
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6 < —40°. The averages listed in columns (7) and (8) were performed over the same 0.5° x 
0.5° pixels for which the average UV intensities were determined. 


The data reduction, point source removal, and the treatment of foreground light sources 
for the FAUST images was described in detail by Sasseen et al. (1995). An additional 
discussion of the impact of the space shuttle environment on the astronomical observations 
and the removal of such effects is provided by Lampton et al. (|1993|) and Chakrabarti et al. 
( p.993|) . With contributions due to direct starlight effectively removed (|Cohen et al. 199^ ), 
the residual UV intensities listed in column (6) of Table 1 should therefore contain mainly 
a sum of DGL, residual airglow, and an essentially isotropic cosmic background. 

Small, additional contributions to the diffuse emissions from the galaxy are likely in 
the form of H 2 fluorescence ( [Martin et al. 19^ ) and of line emission from C IV at 155.Onm, 
Si IV at 139.7nm and O [III] at 166.3nm ( Martin fc Bowyer 1990|) . Martin et al. (1990) 
detected about 75-150 units at low and intermediate latitudes for H 2 fluorescence, while 
Martin & Bowyer found < 20 units for the listed line emissions. These intensities are small 
compared to the variations in the airglow component, and lacking the facility to separate 
them spectroscopically in the FAUST experiment, we will include these emissions in the 
airglow component in subsequent sections and use DGL to refer only to dust-scattered 
radiation. 


We conhrm the presence of a strong DGL component in the data by plotting the 
observed average UV intensities for our 14 target helds against the corresponding HI column 
densities in Figure 1. The DGL component is expected to depend upon the column density 
of dust and upon the illumination conditions of this dust set by the galactic radiation held. 
The column density of dust is correlated with the HI column density ( jBohlin et al. 1975|) . 
Hence, we expect a correlation between the DGL component and HI which Figure 1 indeed 
reveals. 


We also note the interesting fact that the variation in the UV intensity in individual 
helds (Table 1) compared to the variation of the HI column densities for the same helds is 
larger on average (ratio = 3.62 ± 1.37 [photon units/10^® cm“^]) than the same ratio (Fig.l) 
for the total data set (2.50 ± 0.37 [photon units/10^^ cm“^]). This is easily understood if 
one considers the diherence in spatial resolution represented by the two data sets, 30' x 30' 
pixels for the UV intensity measurements vs. the 2° FWHM width of the horn beam of 
the Bell Lab HI survey. We infer that there must be real intensity variations in the UV 
background radiation due to structure in the ISM too small to be resolved by the HI survey. 

The positive intercept at N(HI) = 0 in Figure 1 suggests that the sum of residual 
airglow and extragalactic background is approximately 600 units. It is likely that this 
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Average H I Column Density [10'® cm^®] 


Fig. 1.— A plot of the average diffuse UV background intensity measured in 14 FAUST 
fields against the average HI column densities. Averages are performed over a 4° x 4° area 
centered on the held coordinates. The error bars represent one standard deviation of the 
distribution of values about the respective averages. 


is an underestimate because the coupling of a forward-throwing phase function with the 
concentration of galactic light sources near the galactic plane produces a higher scattered 
light intensity per unit dust column at lower latitudes, i.e., higher N(HI) values, compared 
to higher latitudes, or lower N(HI) values, thus steepening the slope of the correlation and 
lowering the corresponding intercept (WP |1994|) . We can test this by plotting the UV 
intensities against corresponding 100//m IRAS background intensities, as shown in Figure 
2. The IRAS background is correlated both with the dust column density and the density 
of the interstellar radiation held responsible for heating the grains along the line of sight. 
The steepening of the correlation suspected in Figure 1 should therefore be reduced; a 
residual effect could still result from the decrease of radiation density with z-distance from 
the galactic plane. 
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Average IRAS 100 grm [MJy/sr] 

Fig. 2.— The same UV intensities as in Fig. 1 plotted against the diffuse Galactic IRAS 
100;Ltm background intensity. The solid line indicates the linear correlation involving all data 
points; the dashed line shows the linear correlation for the data with the Dorado measurement 
excluded. 


In Figure 2, the intercept at zero IRAS intensity suggests a background of about 835 
units, as suggested above. However, this is rendered somewhat uncertain by three factors: 

(1) the zero-point uncertainties ([Boulanger fc Perault 1988[) of the IRAS lOO/rm detector; 

(2) the likely presence of an as yet undetermined extragalactic background component in the 
100pm intensity; and (3) the exceptional weight of the Dorado measurement on the slope 
of the correlation. Since the high UV intensity at Dorado is consistent with expectations 
from our model (see Figure 5), we suspect that the 100pm intensity of Dorado is in error. 
Eliminating Dorado from the correlation in Figure 2 (dashed line), we arrive at an intercept 
of 550 units, in essential agreement with the intercept in Figure 1. The closeness of the 
correlation between DGL and far-IR background in Fig. 2 is marginally better (r = 0.93, 
Dorado omitted) than the correlation between DGL and HI in Figure 1 (r = 0.79); this 
may be understood through the fact that both DGL and IRAS 100pm intensities depend 













directly on the dust column density, while the DGL-HI correlation involves a possibly 
variable gas-to-dust ratio. We must attribute the real width in the DGL distribution for 
a given IRAS intensity as resulting from variations due to phase function effects in the 
far-ultraviolet, which would leave the IRAS background unaffected, and from deviations of 
the line-of-sight cloud spectrum from the average (see §3.3.4). 

This initial examination of the FAUST data revealed that the measured diffuse 
background intensities consist of a constant component of at least 600 units and a DGL 
component which is correlated with other measures of the galactic dust column density 
and with the galactic radiation density. Provided the hrst component is, on average, the 
same for any random subset of the data, the latter component can be analyzed with an 
appropriate radiative transfer model. 


3. The Radiative Transfer Model 


3.1. General Reqnirements 

A successful model for the DGL must have the following properties: (1) It must contain 
the actual observed illuminating radiation held, (2) the scattering medium should have 
cloud structure, and (3) multiple scattering must be included. These requirements become 
clear when one considers the following points. 

One of the most important characteristics of the DGL at far-UV wavelengths is 
apparent from Figure 1: The DGL intensity for a given N(HI) varies by factors of 
two-to-three with direction in the sky after the approximately constant background has 
been subtracted. While part of this variation may be due to a spatially variable gas-to-dust 
ratio or to dust associated with HII or II 2 rather than HI, most of the observed variation is 
a result of the spatially, severely anisotropic, interstellar radiation held (|Gondhalekar 1990 


Henry et ah 198^ . A hrst requirement of any suitable DGL model for the UV is therefore 
that it must employ the known distribution of the far-UV sources of radiation. A radiation 
held based on the expected UV huxes from known hot stars listed in the SKYMAP star 
catalog ( [Gottlieb 1978| ) has recently been presented by Murthy & Henry (1995). Another 
useful source is the measured UV huxes of some 58,000 stars detected by the TD-1 satellite 
( pondhalekar 19^ , which we have used for this work. 


The interstellar medium responsible for the opacity of interstellar space is distributed 
in the form of structures commonly referred to as clouds. Glouds containing grains with 
single-scattering albedo a < 1 and exposed to an external radiation held will have an albedo 
less than the single-scattering albedo ( [van de Hulst 19871) , and this diherence in albedo will 
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become larger with increasing optical depth of a clond and with increasing phase function 
asymmetry (|Witt fc Gordon 1996| ). In the far-UV the optical depth per cloud is typically 
three times that at visible wavelengths, and previous studies suggest an increase of the 
phase function asymmetry toward shorter wavelengths ( [(Jalzetti et al. 1995 ). Thus, a DGL 
model suitable for the analysis of far-UV data should incorporate a cloud structure for the 
scattering medium if it is to avoid a systematic underestimate of the far-UV albedo. 


Van de Hulst and de Jong ( |1969|) demonstrated that multiple scattering is an 
indispensable feature involved in the transfer of DGL. This applies even more to the 
UV because of the greater optical depths of the galactic dust layer. In the case of the 
present observations, which are restricted to intermediate and high galactic latitudes with 
N(HI) < 8 X 10^° cm“^ (see Table 1), multiple scattering is marginally important. With 
corresponding UV optical depths of order 1.2 or less, the resultant albedo bias can be as 
high as 25% if the entire line-of-sight material is contained in a single cloud and a single 
grain albedo of 0.5 is assumed. We will, therefore, include multiple scattering in our model 


3.2. Model Design 


The principle of the WP radiative transfer model follows the pioneering approach 
of Mattila (1971) by taking into account the discrete cloud structure of the interstellar 
medium. Mattila’s treatment was limited to directions in the galactic plane which placed 
his calculations into the limit of very large optical depths. This meant that the DGL 
intensity was determined mainly by the adopted incident radiation held which was derived 
from photometric measurements of the integrated sky brightness as a function of galactic 
latitude and longitude. 


In the WP model, the approach was generalized to include all galactic latitudes. As 
a source for the (mostly small) optical depths at higher latitudes we adopted the Bell 
Laboratories HI Survey (^tark et al. 19921) , supplemented by the Parkes HI survey (|Gleary 
|et al. 19T9| ; [Heiles fc Gleary 1979|) . We used the conversion from HI column density to 
optical depth at the FAUST effective wavelength of 165nm as given by Sasseen Sz Deharveng 

(EH), 

ri65 = 0.144V^^[10^°cm“^], 


( 1 ) 


which is based on the N(H) - E(B-V) relationship observed by Bohlin et al. (1978) and 
Diplas & Savage (|1994|) , and the average galactic extinction curve of Savage & Mathis 
(1979). 
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3.2.1. Radiation Fields 

In order to extend the radiative transfer model to directions with |6| > 0°, we need 
to express the variation of the radiation held with z-distance from the plane. We start by 
constructing the radiation held in the galactic plane (z = 0) utilizing the huxes of about 
58,000 stars measured by TD-1 at A = 156.5nm (FWHM = 330A) which were kindly 
provided to us by J. Murthy. Maps of this radiation held may be found in the recent work 
of Murthy Sz Henry (1995) and in the Atlas of the Ultraviolet Sky by Henry et ah (1988). 

The TD-1 catalog is thought to be complete down to huxes of 1.0 x 10“^^ erg cm“^ s“^ 

at 156nm. The contribution by stars fainter than this limit to the integrated stellar 
huxes appears to be negligible at this wavelength, certainly much less than 10% of the total 
(Gondhalekar 1990). However, the DGL observed by TD-1 ([Morgan et ah 1978| ) is part of 
the illuminating radiation held seen by each cloud and needs to be added to the integrated 
stellar hux. We therefore computed the expected DGL radiation held for a dust albedo 
a = 0.5 and a strongly forward directed phase function and added it to the integrated 
stellar huxes. This increased the total background intensities by typically 25%. 

For the purpose of our calculation, we represented this radiation held (at z = 0) as 
an 800-element, two-dimensional array, and expressed the average intensity in each array 
element in photon units. The spacing of the array elements and their sizes are non-uniform. 
At low galactic latitudes, and other directions where large gradients in the radiation held 
intensity occur, elements are spaced more closely, while at high galactic latitudes large 
elements are chosen. 

The transformation of the central-plane radiation held into one seen by a cloud 
located at distance zc % 0 from the galactic plane was approximated following a method 
described by Mattila ( |198U4 [SD. For the purpose of this transformation only, we assumed 
that the principal source of the A = 156.5nm radiation held is galactic OB stars whose 
volume emissivity is distributed in a plane-parallel disk with an exponential z-distribution 
characterized by a vertical scale height /3 = 60 pc. We assume this because about 90% of 
the total integrated stellar radiation at 156nm arises near the galactic plane dehned by 
—45° < b < 45° (Gondhalekar 1990). Also, this assumption appears justihed because the 
number of suspected subdwarfs in the TD-1 catalog, though greater than anticipated, was 
only of order 10^ ( jGarnochan fc Wilson 19831) . The galactic distribution of these stars does 
not diher drastically from that of the overall stellar radiation held and their contribution to 
the radiation held was included in any case. By comparing the expected 5-distribution of the 
integrated intensity seen by an observer at z = zc with the symmetric distribution seen at z 
= 0, we derived transformation factors as a function of zc and b. This was done by dividing 
the skewed radiation held distribution at z = zc by the symmetric distribution at z = 0. To 
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obtain an actnal radiation field applicable to z = z^, we then mnltiplied the TD-1 radiation 
field with the appropriate ^-dependent transformation factors. As a resnlt, radiation fields 
at increasing \zc\ exhibit a growing asymmetry in the radiation emanating from the two 
galactic hemispheres with b > 0° and 6 < 0°, and the near-eqnatorial intensity peak shifts 
toward latitudes in the opposite hemisphere while maintaining the basic non-isotropy and 
longitudinal variability of the original TD-1 field. We constructed a total of 17 radiation 
fields, for values of zq = 0, ±25pc, ±50pc, ±75pc, ilOOpc, ±150pc, ±200pc, ±250pc, and 
±300pc. 


3.2.2. Cloud Spectrum 


Interstellar clouds exist with a wide range of optical depths (e.g. Crovisier 1981). In 
order to reflect these conditions we adopted a cloud spectrum of three types of spherical 
clouds: uniform clouds of optical depth diameter (at 156.5nm) of 0.5, uniform clouds of 
optical depth diameter of 2.5, and centrally condensed clouds (p oc r“^ for r/R > 0.05; 
p = constant for r/R < 0.05) of optical depth diameter of 10. The probability ratio of 
encountering such clouds along any line of sight was taken to be 12:3:1. Assuming arbitrary 
impact parameters for a penetrating line-of-sight, such a cloud spectrum provides, on 
average, a column density distribution very similar to that derived by Crovisier ( p.981|) 
from the statistical properties of interstellar HI which were based on 21-cm line absorption 
surveys. The average optical depth per cloud for our adopted spectrum was found to be 
0.593 by subjecting our spectrum of three clouds, in a frequency ratio of 12:3:1, to random 
penetration by a line of sight. The Crovisier study seemed to be particularly relevant in the 
current context: (1) its database consists mainly of observations out of the galactic plane, 
as is the case with the FAUST data; (2) our source for the column density of dust is the 
correlation between E(B-V) and HI ( [Bohlin et al. 1978[) ; (3) the range of total HI column 
densities covered by the Crovisier study coincides with the range encountered in the FAUST 
fields in this study. 


The probability of encountering a cloud of any type at a distance z from the galactic 
plane was assumed to be given by an exponential distribution with a scale height of 110 pc, 
i.e. was assumed to be proportional to exp (-z/110 pc) (Crovisier 1981; [Lockman 1984|). 
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3.3. Model Operation 

The computation of the expected DGL surface brightness in a given direction with our 
model involves several Monte Carlo processes: the multiple scattering radiative transfer 
within each cloud of our cloud spectrum, the determination of the number of clouds in a 
given direction, the determination of the type and location of the cloud along the line of 
sight, the determination of the individual surface brightness prohle of each cloud along the 
line of sight, the determination of the penetration point of each cloud, the determination 
of the optical depth along the chosen line of sight through each cloud, and hnally the 
integration of intensities of the clouds along the line of sight. 


3.3.1. The Radiative Transfer in Individual Clouds 

The transfer solution for a spherical cloud of given optical depth diameter, dust 
scattering parameters given by albedo a and phase function asymmetry g, and density 
prohle was found with the method outlined in Witt & Stephens (1974) and in FitzGerald 
et al. ( p.97ti|) . The method traces the transfer of individual photons in reverse and produces 
weight matrices which express the probability that photons hitting the cloud from any 
direction will exit the cloud at a given projected radius in the direction of the observer. 
For the non-isotropic scattering phase function, the Henyey-Greenstein ( p.941|) function was 
chosen. Our approach has the advantage that the radiative transfer through individual 
clouds needs to be done only once for a given set of parameters. The specihc surface 
brightness prohle exhibited by a cloud at a given zq can be obtained by multiplying the 
cloud’s probability matrix with the appropriate radiation held matrix corresponding to the 
value of zc. 


3.3.2. The Number, Type, and Location of Clouds along the Line of Sight 


The average number of clouds along a given line of sight with total column density 
N(HI) is found by dividing the corresponding optical depth derived from Eq. (1) by 
the average optical depth per cloud < tq > = 0.593 (§3.2.2). For a typical line of 
sight in the galactic plane (Ay = 1.9 mag/kpc; [Spitzer 1973|), with an average ratio of 


N(HI)/E(B-V) = 4.8 X 10^^ atom cm ^ mag ^ ( [Bohlin et al. 1978|) , we hnd 7.1 clouds/kpc. 


in close agreement with estimates discussed by Spitzer (1978). For typical high-|6| lines of 
sight, we generally hnd 2, 1, or no clouds. 


The actual number of clouds, n, for a given Monte Garlo simulation is then selected 
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by a Monte Carlo process from a Poisson distribution whose average is the average 
number of clouds just determined. The types of clouds (optical thickness diameter, density 
distribution) are found by random-number choices from within our cloud spectrum based 
upon their relative frequencies. The z-distances of these clouds are then derived from yet 
another series of Monte Carlo processes applied to the assumed exponential distribution 
of cloud distances from the Galactic plane. This sampling of cloud distribution by 
number, types, and location along the line of sight was performed 300 times for each of the 
approximately 60 pixel positions within each FAUST field for a total of 1.8 x 10^ integrations 
per field. Thus, the statistical properties of the interstellar medium are interwoven into the 
prediction of the DGL intensity. We therefore predict an actual distribution of expected 
DGL intensities for each pixel, whose FWHM provides a measure of the DGL uncertainty 
resulting from possible deviations from average conditions. Examples of such distributions 
are shown in Figure 3 and will be discussed below in §3.3.4. 


3.3.3. Surface Brightness Profiles and Total DGL Intensities 


The probability weight matrix of each of the randomly selected and positioned clouds 
(§3.3.2) must be multiplied by the appropriate radiation held in order to predict their 
radial surface brightness prohles. The most abundant variety of clouds, the optically thin 
ones, exhibit a brightness prohle with a maximum at the projected center; the moderately 
optically thick clouds present a fairly hat surface brightness prohle with limb darkening, 
while the highly optically thick clouds display a dark central region surrounded by a bright 
rim ([Witt fc Stephens 1974]) . The surface brightness value added by each cloud along the 
line of sight is determined by the selection of a random penetration point on the geometric 
cross section of a cloud. In order to correct the surface brightness of still more distant 
clouds seen in this direction, we also determine the optical depth through the entire cloud 
at the penetration point. If the intensities of the clouds along a given direction are denoted 
by Ii, U, I3, ..., the total integrated intensity for this line of sight ([Mattila 1971| ) is 

DGL = h + he-^^ + + ... 


I 


( 2 ) 


3 . 3 . 4 . Dependence of Results on Model Parameters 

In its attempt to represent the DGL radiative transfer in as realistic a manner 
as possible, the model described above involves numerous parameters. While most of 
these, such as the radiation field distribution and the optical depth distribution, are well 
constrained by available observations and may therefore be considered fixed, it is the nature 





14 


of the adopted cloud spectrum that affects the model predictions most profoundly. As 
indicated in §3.1, the greater the optical thickness of an individual cloud, the lower its 
resultant effective albedo for a given radiation field and single grain albedo. In addition, 
in directions with low column densities, a cloud spectrum incorporating clouds of higher 
optical depth naturally leads to instances of ’’holes”, i.e. directions with no clouds and 
thus, zero DGL intensity, because the beam-averaged optical depth must still be that 
constrained by the HI column density. If, on the other hand, the cloud spectrum consists 
entirely of very-low-r clouds, there will always by numerous clouds along lines of sight with 
small overall column density, assuring a more uniform sky coverage with more efficiently 
scattering clouds. 

We illustrate these effects in Figure 3 by showing the results of 100 samplings with 300 
photons each along a single direction {i = 315°, 6 = -|-35°, N(HI) = 5.15 x 10^°cm“^) for 
three different cloud spectra. The hrst consists of optically thin clouds of type 1 only, with 
an optical depth diameter of 0.5 at 156.5nm wavelength. The average optical depth for 
cloud penetration is 0.33. The predicted DGL intensity is 1871 ± 59 units. The second is 
our adopted case of a cloud spectrum with ratios 12:3:1 divided among clouds of type 1, 2, 
and 3, respectively. The average optical depth per cloud is now 0.59 and the predicted DGL 
intensity is 1329 ± 52 units, some 30% lower than case 1. The third case displayed is for 
a cloud spectrum 1:1:1, which now emphasizes optically thick clouds. The average optical 
depth per cloud is 0.83 and the predicted DGL, as expected, is still lower at 1101 ± 52 
units. All calculations assumed the same scattering characteristics, namely a = 0.45 and 
g = 0.675. 


For the analysis of the FAUST data we have adopted a cloud spectrum which is 
dominated by low-r diffuse clouds, but also contains a small number of clouds of greater 
optical thickness in a ratio of 12:3:1 in terms of their average relative frequency along a 
given line of sight. This spectrum reflects the demonstrated presence of numerous small 
molecular clouds at high galactic latitudes ([IVfagnani et al. 1986|; |lVfagnani et al. 1996|; |Stark 


1995 ; see Magnani 1994 for a recent review). Our adopted spectrum is also consistent with 


the observed distribution of HI column densities based on the 21-cm line absorption survey 
by Grovisier ( [1981|) . It should be clear from Fig. 3 that the albedo derived from the FAUST 
data using our adopted cloud spectrum will be higher by about 30% than it would be, had 
we assumed low-r clouds entirely. The dependence of the intensity of the predicted DGL 
upon the chosen cloud spectrum must also be kept in mind when we compare predicted 
DGL intensities with those actually observed (see Fig. 6). 
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0 1000 2000 3000 

Average Intensity [photon units] 

Fig. 3.— Model prediction of the intensity of UV DGL plus UV background derived from 100 
samplings with 300 photons each for three different cloud spectra. All calculations are for 
the identical direction [i = 315°, b = +35°) and for identical grain properties (a = 0.45, 
g = 0.675). The hydrogen column density is N(HI) = 5.15 x lO^'^cm”^. 

4. Derivation of Scattering Properties 

Our aim in this section is two-fold: (1) We want to separate the galactic component of 
the diffuse UV intensities measured by FAUST (Table 1) from the non-galactic components. 
Contributions to the latter are terrestrial airglow and extragalactic background radiation. 
(2) We want to analyze the galactic component using the DGL model described in the 
previous section. The aim is to derive the average dust albedo a and the average phase 
function asymmetry parameter g which characterize the galactic dust in the diffuse ISM 
at intermediate and high galactic latitudes. Toward this end we calculated the expected 
average UV intensity for each of the 14 FAUST helds, each based on approximately 1.8 x 10"^ 
integrations, for a wide grid of closely spaced values of a and g. 










16 


4.1. The Fitting Criteria 


Which set of scattering parameters for dust grains will lead to the best representation 
of our data? To decide this question, we assumed the following, (a) All FAUST lines of 
sight contain dust with identical average properties. Since for our observations |6| > 21°, 
most of the clouds studied are well outside the galactic plane and are representative of 
the low-density diffuse interstellar medium, (b) The measured surface brightness of each 
FAUST held contains an unknown component with a surface brightness > 0 units due to 
residual airglow and a component of extragalactic origin. We assume the latter to have an 
intensity of 300 units (WP |1994|) outside the Galaxy. This component is assumed to be 
isotropic but to be attenuated for each line of sight by the galactic dust column derived 
from the N(HI) values. We included this extragalactic component both into the illuminating 
radiation held (essentially negligible) and into the predicted intensity for each FAUST held, 
(c) On average, the FAUST helds at lower galactic latitudes are subject to the same airglow 
contamination as the FAUST helds at higher galactic latitudes. This does not require that 
each held have an identical airglow component, only that any random group of six or seven 
helds have similar airglow averages. Thus, the galactic latitude dependence seen in our 
measured intensities is due entirely to DGL and the extinction-modulated extragalactic 
background with the former generally increasing with decreasing latitude and the latter 
decreasing with the increasing extinction at lower latitudes. Given these assumptions, the 
best-ht model must satisfy these requirements: 


1. Slope Griterion: The model predictions of the DGL values for the FAUST helds, 
including the partially attenuated extragalactic radiation, must exhibit the same 
slopes as the observations when plotted against the corresponding values of N(HI), 
cosec|6|, and the IRAS lOO/rm intensities for the FAUST helds. 

2. Positive Airglow Griterion: The average predicted intensities of DGL and extragalactic 
light for individual FAUST helds should in no case fall above the measured intensities 
since this would imply negative airglow. 

3. Equal Airglow Griterion: A considerable range of dust scattering parameters {a,g) 
will be able to satisfy the hrst two criteria. However, for most sets of values of {a,g), 
the value of the positive airglow derived from application of criteria (1) and (2) is 
diherent when observations and model predictions are plotted against N(HI) and when 
plotted against cosec|f)|. The equal airglow criterion demands that the value of the 
airglow deduced from the comparison of intercepts is the same when observations and 
model predictions are plotted against N(HI) and cosec|6| and evaluated at N(HI) = 0 
and cosec|6| = 0. 
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A much more limited set of values (a, g) will satisfy this third criterion in addition 
to satisfying the first two. We do not include the intercept in the plot against the IRAS 
100/xm diffuse intensity in this criterion because calibration uncertainties in the zero point 
of the IRAS intensity scale and the presence of some extragalactic component in the IRAS 
data would render the outcome uncertain. 


4.2. Models 

In Figure 4, we show all combinations of a and g for which complete models were 
calculated. Different symbols indicate the slope found when predicted intensities are plotted 
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Fig. 4.— Complete models were calculated for the grain scattering properties indicated by 
the symbol locations in the (a, g) plane. Different symbols indicate the slopes in graphs of 
the predicted UV intensities vs. N(HI). The encircled symbols correspond to models which 
predict intensities higher than observed. 

against the N(HI) values for the corresponding FAUST fields. The observed slope (see Fig. 
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1) is 2.5 [photon units/10^®cm“^]. It is clear from Fig. 4 that the htting criterion (1) is met 
within a very narrow range of a-values (0.4 < a < 0.5) but for a large range of gf-values. 
The positive airglow criterion eliminates all models identihed by encircled symbols, as well 
as all those with still larger a-values which were left uncomputed. 

As an example, a model with a = 0.45 and g = 0.675 is shown in Fig. 5 in comparison 
with the FAUST data plotted against the N(HI) values of the corresponding FAUST helds. 
The model points, denoted by open circles with vertical error bars, follow the observed 



Fig. 5.— Model predictions (open circles) are plotted together with corresponding 
observations against average N(HI) values for FAUST helds. The vertical model error bars 
denote one standard deviation in the distribution of predicted intensities for a 4° x 4° held. 

intensities closely, except they are displaced by a constant vertical ohset attributed to 
the airglow component. For example, the model successfully reproduces the signihcant 
diherences in intensity observed for groups of FAUST helds for which the HI column density 
is almost constant. Specihc cases are the groups consisting of Hydra 21c, Centaurus and 
Hydra 21b, and Virgo PI, Dorado, and Virgo P2. This gives us conhdence both in the 
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model’s predictive power and in our assumption that the residual airglow contribution is 
roughly constant. Note that the positive airglow criterion (2) is also met in the model 
shown in Fig. 5. 

To further illustrate the close match between the model predictions with the individual 
observations, we added an average airglow intensity of 530 units to the model values and 
overplotted these on the observations in Figure 6. In almost all cases, the observations are 



Fig. 6.— The same data and model predictions as in Figure 5, except that a constant airglow 
intensity of 530 units has been added to the model points. Error bars have been omitted for 
clarity; they would be identical to those shown in Figure 5. Note that the observed points 
for NGC 6752 and Hydra 21a essentially coincide. The corresponding model points are those 
immediately above and below the observed points. 


reproduced within better than one standard deviation. 
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4.3. The Equal Airglow Criterion Applied 

Thus far we have found combinations of a and g which reproduce the slope of the 
observed data in Figure 5. 

As mentioned, we interpret the vertical offset between the least-squares slopes htted 
to the observed points and the predicted points in Fig. 5 as resulting from airglow. The 
airglow intensity derived from Fig. 5 is about 530 units, on average. When the same data 
and predictions are plotted against cosec|6| of the positions of the FAUST helds, we arrive 
at Fig. 7. Again, the observed and the model slopes are nearly identical and the vertical 



Average Cosee lb 


Fig. 7.— The same data and models as in Figure 5 plotted against cosec \b\ of the FAUST 
helds. The “best ht” solution shown here requires that the vertical offset between observed 
slope and model slope be the same as in Fig. 5. 


offset between the two slopes is approximately 530 units, i.e., the inferred airglow intensity 
is the same as derived from Fig. 5. If the gf-values of the models are increased beyond 0.675, 
the cosec 15| plot yields higher airglow values than the N(HI) plot as well as divergent slope 
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values; if g is reduced below 0.675, the reverse occurs. This is shown quantitatively in Fig. 
8 where lines of constant airglow as derived from the two types of plots are shown in the 
relevant area of the (a, g) plane. Allowable valnes for a and g must be located in the part 
of the plane where the two comparisons yield the same airglow values. 



Fig. 8.— Contours of constant airglow dednced from comparing models and observations 
when plotted against N(HI) and against cosec|&|. This graph snggests an average airglow 
near 510 photon units as being associated with the “best £t” solution. 


The straight forward application of the three fitting criteria constrains the best solution 
to a narrow region surronnding the point a = 0.45 and g = 0.675 in the {a^g) plane. Using 
this solntion, we have derived the likely airglow intensities for our 14 FAUST fields and 
have listed them in Table 2. 





Table 2. Deduced Airglow Values^ 


FAUST Field 

Deduced Airglow 


(Photon units) 


Dorado 

582 

± 

233 

NGPKM 

400 

± 

70 

M 87 

883 

± 

93 

Virgo P2 

470 

± 

79 

Virgo PI 

559 

± 

142 

Centaurus 

657 

± 

368 

M 83 

449 

± 

368 

NGC 6752 

457 

± 

344 

Hydra 20a 

307 

± 

298 

Hydra 20b 

307 

± 

220 

Hydra 20c 

390 

± 

239 

Hydra 21a 

863 

± 

263 

Hydra 21b 

717 

± 

338 

Hydra 21c 

494 

± 

391 


^Based on the assumption that the 
extragalactic isotropic background before 
galactic extinction is 300 photon units. 
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4.4. Uncertainties 


It was shown in Fig. 4 that the albedo value of the ’’best £t” solution is sensitively 
dependent upon the slope of the measured FAUST diffuse intensity plotted against N(HI). 
Because of the close correlation between N(HI) and cosec|6| and between N(HI) and the 
IRAS lOOpm diffuse galactic background ( [Boulanger fc Perault 198^; [Boulanger et ah 


19961) , the same would be the case, if plots are made against either of these two other 


quantities. 

There are two major sources of uncertainty in these observed slopes: (1) The residual 
airglow is variable with time and direction, and thus, a randomly variable component 
is added to the true DGL values. The values in Table 2 are representative of the range 
of variability to be expected for the airglow component. Only if our DGL model were 
perfectly accurate would Table 2 give the actual airglow values, and we do not assume 
this. The presence of a random airglow component having the same distribution as Table 
2 among our data renders the slope of the linear regression in the plot against N(III) 
uncertain by ± 0.34 [photon units/10^®cm“^]. (2) The absolute calibration of the FAUST 
instrument is uncertain by about 15% (Powyer et ah 1993[) . While this does not affect 
the uncertainty of the linear regression, it affects each data point by the same percentage, 
and thus both the intercept and the slope are uncertain. We hnd the impact of the 
systematic calibration uncertainty upon the slope to be ± 0.37 [photon units/10^®cm“^]. 
Adding these uncertainties in quadrature, we hnd for the observed slope 2.5 ± 0.5 [photon 
units/10^®cm“^]. 

An additional uncertainly arises from the fact that two independent model calculations 
with identical parameters but different seeds for the random number generator produce 
model slopes which typically vary by ± 0.10 [photon units/10^®cm“^], which must be added 
to the uncertainty above when observations and models are compared. 


4.5. Scattering Properties 

We again plot in Fig. 9 the (a, g) plane as in Fig. 4 in which the region of possible (a, g) 
values is now delineated by the application of the three htting criteria and the uncertainties 
in the slope discussed above. The outer boundary is set by the slope uncertainty in the 
plot of FAUST intensities against N(HI). The cross-hatched area is considered unlikely for 
possible solutions since the model, with no airglow, overpredicts the observed intensity. The 
equal airglow criterion and the slope criterion lead to a most probable value for the albedo, 
a = 0.45 ± 0.05, and for the phase function asymmetry, g = 0.68 ± 0.10. 
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Fig. 9.— The derived scattering properties in the {a,g) plane as constrained by the three 
htting criteria (§4.1). The slopes referred to are those of the predicted UV-DGL-N(HI) 
relationships. 


5. Discussion 


5.1. Scattering Properties 


Ultimately, we would like to know whether the population of dust grains is different in 
different astrophysical environments within the Galaxy. Standard dust grain models are 
based mainly on the observation of the average galactic extinction curve (Savage & Mathis 
1979) and observed depletions of heavy elements from the gas phase (|Soha et ah 1994| ; 
[Snow fc Witt 19"^ |1996|) , and are therefore representing mostly the diffuse interstellar 
medium near the Galactic plane. However, differences in extinction characteristics between 
different lines of sight (|Witt et al. 198^ ; [Fitzpatrick &: Massa 199U| [Mathis fc Gardelli 


1992[; [Jenniskens fc Greenberg 1993[) as well as differences in depletion patterns ([SembacE 


fc Savage 1996[) are well documented, and these are reflected in differences in the size 
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distributions of grains derived for different environments (e.g. |Kim et al. 1994|) . 

Scattering characteristics, like extinction properties, are also dependent on grain 
composition and size. There is little reason to suspect that dust grains in different 
environments, e.g. reflection nebulae, HII regions, or the low-density diffuse ISM, are 
similar with respect to their scattering properties. This is especially so in light of the 
considerable observational evidence that the extinction properties of the dust exhibit 
substantial variations in the UV, both in magnitude and in wavelength dependence. 

In Table 3 we have compiled a number of values of (a, g) which have been derived 
recently from the observation and analysis of scattered UV light in reflection nebulae and 
star forming regions, all regions characterized by densities in the range 10^ cm’ -3 < n(H) 

cm“^. The denser reflection nebulae appear to share albedo values in excess of 0.6 in 
the 130-170nm spectral range, while the present study yields a value dehnitely lower than 
that. This difference is most easily understood if we assume that the average grain size is 
larger in reflection nebulae, or, what would amount to the same, that the number of very 
small grains (size A) is reduced compared to the number of larger (size ~ A) grains. This 
is indeed supported by the fact that in reflection nebulae Ry = Ay/E(B-V) is usually found 
to be larger than the ’’normal” galactic value for the diffuse ISM of Ry = 3.1, which is true 
also for the regions compared in Table 3. In addition, most of the dust observed in this 
study is seen at relatively high galactic latitudes, placing it potentially at large z-distances. 
With our assumed z-distribution, 50% of the dust is located at |z| > 75 pc. Sembach and 
Savage (1996) have shown that dust grains in such environments are subject to erosive 
processes which return a portion of their solid substance to the gas phase. Kim and Martin 
( p.996|) have shown that size distributions of interstellar grains characterized by Ry >3.1 
have higher albedos than the standard Ry = 3.1 grains throughout the entire near-IR to 
UV spectral range. There are still discrepancies, however, between the observationally 
derived albedo values for reflection nebulae and the values predicted by Kim & Martin. 
This may simply reflect the fact that their assumed graphite/silicate composition for the 
grains does not correspond entirely to reality. Nevertheless, we have for the hrst time a 
relatively clear indication that the dust albedo in the UV in the diffuse ISM of the Galaxy 
is lower than the albedo in denser regions, about 0.45 vs. 0.65, and that this difference is 
related to the typical grain sizes found in these environments, with smaller grains in the 
diffuse ISM yielding a lower albedo. 

Additional support for these conclusions is found in the observations and analysis of 
scattered far-UV radiation seen at high galactic latitudes by Onaka and Kodaira (|1991|) . 
While more limited in spatial coverage, their data set allowed them to conclude a > 0.32 
and g > 0.5, in agreement with our results. The suggestion of smaller dust grains, on 
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average, at higher |z|- distances also finds support in the finding by Kiszkurno & Lequeux 
( p.987| ) of a general steepening of the far-UV extinction curve with increasing |z|-distance 
from the galactic plane. This characteristic is attributed to a systematic shift toward 
smaller grain sizes by these authors. 

When comparing the phase function asymmetry values in Table 3, no clear difference 
can be established, and, in fact, may not exist. Within existing uncertainties, we conclude 
that the phase function asymmetry of scattering by dust in the diffuse ISM and in dense 
nebulae is identical. The phase function asymmetry is determined solely by those grains 
which contribute to the scattering. If, for example, the diffuse ISM contains a separate 
population of very small grains (-C A) which would provide absorption, the albedo would be 
lowered while the phase function would still entirely be determined by the larger scattering 
grains and would remain unchanged. In addition, a very small change in g represents a 
significant change in the phase function when g 0.7, and as long as any determination 
of g is uncertain to ± 0.1, any real differences equal to small A 5 ^’s will not likely be 
detectable. The only situation which allows more precise differences in g to be derived is 
the determination of relative ^f-values as a function of wavelength in a given object (e.g. 
Calzetti et ah 19951) . 


5.2. The DGL-N(HI) Correlation 

The positive correlation between the far-UV background intensity and the column 
density of neutral hydrogen is generally viewed as evidence for the Galactic origin of a 
major part of the UV background. ( [Lequeux 199CI| ; Bowyer 1991). However, a number of 
separate investigations (Bowyer 1991, Table 1) have produced a wide range of values for 
the slope of the DGL-N(HI) correlation, often mutually exclusive. Bowyer ( |1991|) tabulates 
values ranging from 0.3 to 1.8 [photon units/10^® HI cm“^] from eleven different studies. To 
this, we now add our present result of 2.5 ± 0.5 [photon units/10^® HI cm“^], larger than 
any found previously. 

The significance of these apparent discrepancies is easily understood. The scattered 
light intensity per [10^® HI cm“^] is determined by the scattering cross section per H-atom, 
the incident radiation field, and the phase function asymmetry g. A relatively large value 
of g, as we have deduced in this investigation, provides a strong coupling between the 
anisotropy of the illuminating radiation field and the anisotropy of the scattered radiation 
field. The FAUST observations were conducted in a region of the UV sky where the gradient 
in illuminating intensity between near-equatorial galactic latitudes and polar latitudes is 
one of the strongest in the entire sky. The fact that this strong gradient is indeed reflected 
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in a larger value for the slope in the DGL-N(HI) correlation is by itself an independent 
indication that the scattering phase function must be relatively large. Correlations between 
measured UV background intensities with N(HI) are therefore comparable only, if the data 
have been obtained from identical regions of the sky. 

We can demonstrate that our value for the slope of the DGL-N(HI) correlation is well 
within the range expected for different regions of the sky with the following experiment. For 
this test, we adopted scattering properties a = 0.45, g = 0.6 and the average N(HI) values 
of the FAUST helds. We predicted the observable DGL, including 300 units of extragalactic 
background, for the actual FAUST helds and found a slope for the DGL-N(HI) correlation of 
2.58 [photon units/10^® HI cm“^] (Fig. 10). Retaining the galactic latitudes of the FAUST 
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Fig. 10.— Predicted correlation between UV DGL and N(HI) for two different Galactic 
regions, assuming a = 0.45, g = 0.60; the actual positions of the FAUST helds (hlled circles; 
slope = 2.58 [photon units/10^®cm“^]) and hypothetical positions with the Galactic latitudes 
and N(HI) values identical to the FAUST positions, but moved to I = 145° (open circles; 
slope = 0.72 [photon units/10^®cm“^). 




fields as well as their N(HI) values, we calculated the DGL (plus 300 units extragalactic 
background) under the assumption that their galactic longitude was £ = 145°. This is a 
region where the illuminating radiation held has a much weaker gradient with latitude. The 
resulting slope in the DGL-N(HI) correlation was now 0.72 [photon units/10^® HI cm“^], 
near the lower end of the range tabulated by Bowyer (|1991D . This clearly identihes the 
anisotropy of the illuminating radiation held as the dominant cause of the anisotropy of 
the DGL in the presence of a strongly forward-directed scattering phase function. The high 
value of the slope of the DGL/IRAS lOOpm correlation found by Sasseen & Deharveng 
( p.996|) for these FAUST data is also consistent with the picture of anisotropic illumination 
coupled with a high value of g. 


5.3. The Extragalactic Background 


The present study did not permit a separate determination of the airglow and the 
partially attenuated extragalactic background components, except that their sum should 
be about 700 photon units averaged over the 14 FAUST helds. The reason for this is 
two-fold: We lacked a spectroscopic capability which would have been able to separate the 
known line emissions of the airglow from the quasi-continuous extragalactic background 
spectrum; and the variations in airglow intensity among different FAUST helds were large 
in comparison with the weak Galactic latitude dependence in the extragalactic background 
imposed by Galactic foreground extinction. We assumed an extragalactic component of 
300 units in agreement with the determination by WP (|1994| ) and the summary of earlier 
results by Paresce ( p.990|) . This produced the deduced airglow intensities listed in Table 2. 
After correction for Galactic foreground extinction, the extragalactic component is reduced 
to about 230 units for the high-latitude helds and to about 100 units for the FAUST helds 
with the highest N(HI) values leading to an average for our sample of 162 ± 46 photon 
units. 


A discussion of the origin of the extragalactic component and the potential contributors 
is beyond the scope of the present paper. However, we must examine whether the 
assumption of an extragalactic component, anticorrelated with N(HI) through extinction, 
has a signihcant ehect upon the hnal scattering characteristics derived from the DGL 
analysis. Repeating our analysis with the assumption that an extragalactic component 
in the UV is completely absent, we hnd that the excess intensity must then be mostly 
attributed to a higher airglow. The resultant albedo and phase function asymmetry 
are essentially unchanged, 0.43 instead of 0.45, and 0.65 instead of 0.675, respectively. 
Thus, should future studies reveal our assumption of 300 units for the extragalactic UV 
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background as inaccurate, this would not affect the validity of the present results. 

On the other hand, an extragalactic component of more than 300 units is effectively 
excluded by the lower limit to the total far-UV background, including extragalactic, 
measured by Onaka and Kodaira (|1991|) who directly recorded intensities as low as 300 units 
at high galactic latitudes, effectively using the same bandpass. It is important to state, 
however, that it is not necessary to know the intensities of the airglow and the extragalactic 
background separately in order to determine the dust scattering properties of galactic dust, 
provided two conditions are met: (1) the observations must cover a sufficiently large range 
of galactic latitudes to allow the dependence of the measured intensity upon the HI column 
density to be revealed, and (2) the airglow component must be relatively constant during 
the observation. 


6. Conclusions 

We have analyzed the 140-180nm UV background intensity in 14 FAUST fields at 
intermediate and high Galactic latitudes after contributions from discrete sources such as 
stars and galaxies had been removed. We conclude the following: 

(1) We confirm that a major component of the measured UV background must be of 
Galactic origin based on its strong correlation with N(HI), with cosec |f)|, and with the IRAS 
lOO/Ltm background intensity. We interpret this component as DGL, produced by scattering 
of stellar photons by Galactic dust. 

(2) The slope in the FAUST DGL-N(HI) correlation of 2.5 ± 0.5 [photon units/10^® HI 
cm“^] is found to be steeper than that seen in earlier investigations in other Galactic 
regions. We explain this discrepancy as owing to the existence of a very strong Galactic 
latitude gradient in the illuminating radiation field in the FAUST region. 

(3) We found our radiative transfer model capable of matching the detailed variation 
of the DGL intensity from one line-of-sight to another. The most important characteristic 
of the model responsible for this is the incorporation of the detailed anisotropy of the 
interstellar radiation field in the UV and its variation with z-distance from the Galactic 
plane. 

(4) The radiative transfer analysis of the DGL component with our model resulted in 
a set of well-constrained scattering parameters for dust in the diffuse ISM at intermediate 
and high Galactic latitudes: albedo a = 0.45 ± 0.05 and phase function asymmetry 

g = 0.68 ± 0.10. 
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(5) The dust albedo is some 50% lower than that commonly derived for dense reflection 
nebulae and star forming regions. We interpret this difference as arising from a difference 
in size distributions with the grains in the diffuse ISM being smaller on average. 

(6) A contribution of about 700 ± 200 photon units appears to be uncorrelated with 
Galactic lines-of-sight. We interpret this component as due to a sum of residual airglow 
(530 ±190 photon units) and isotropic extragalactic background radiation (160 ± 50 
photon units, after correction for Galactic extinction). 

We thank Jens Petersohn for his invaluable assistance with the WP radiative transfer 
model during the early stages of this work. The referee, Dr. Ken Sembach, provided a 
number of thoughtful and constructive comments which helped to improve the presentation 
of this paper, and for which we are grateful. ANW and BGF acknowledge material support 
from NASA LTSA Grant NAGW-3168 to The University of Toledo. 



Table 3. 


Region 

A (nm) 

a 

9 

Ry 

Ref 

NGC 7023 

144,152 

> 0.65 ±0.09 

0.75 ± 0.05 


1,5 

NGC 7023 

130 

0.62 ± 0.09 

0.75 ± 0.05 


2,5 

Sco OB2 Assn 

177 

0.64 ± 0.09 

nndet. 

3.5 

3 

IG 435 

156 

0.75 ± 0.05 

0.7 ± 0.1 

5.3 

4,6 

Diffnse ISM 

160 

0.45 ± 0.05 

0.68 ± 0.1 

3.1 

this stndy 
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